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Abstract — Advantages  and  limitations  of  washout  filters 
in  feedback  control  of  both  continuous-time  and  discrete¬ 
time  systems  are  discussed  and  generalizations  that  alleviate 
the  limitations  are  presented.  Some  previously  unpublished 
results  in  the  Ph.D.  dissertation  of  one  of  the  authors  (Lee, 
1991)  are  presented  in  the  context  of  their  relation  to  the 
generalized  results  and  to  recent  publications  on  delayed 
feedback  control.  We  show  that  delayed  feedback  control  (for 
discrete  time  systems)  extensively  used  in  control  of  chaos  is 
a  special  case  of  washout  filter-aided  feedback.  Moreover,  the 
limitations  of  delayed  feedback  control  can  be  overcome  by 
the  use  of  washout  filter-aided  feedback,  which  gives  rise  to 
the  possibility  of  stabilizing  a  much  larger  class  of  systems. 

I.  Introduction 

It  is  a  common  practice  in  the  analysis  of  nonlinear 
systems  and  in  feedback  control  design  to  assume  that  the 
equilibrium  point  (or  the  operating  point)  of  the  system  is 
accurately  known  or  does  not  change  over  the  operating 
regime.  However,  models  of  physical  dynamical  systems 
are  in  general  uncertain.  Therefore,  static  feedback  control 
is  ineffective  in  addressing  problems  where  the  operating 
point  is  not  accurately  known  or  there  is  parameter  drift. 

Consider  the  nonlinear  system  described  by 

x  =  f(x,u)  (continuous-time)  (1) 

or 

a;(fc+l)  =  f(x(k),u(k))  (discrete-time)  (2) 

where  /(•,  •)  is  uncertain,  u  is  the  scalar  input  and  x  £  3?”  is 
the  state  vector.  Due  to  the  uncertainty  in  /,  the  equilibrium 
points  (if  any)  of  the  system  (1)  and  the  fixed  points 
(if  any)  of  (2)  are  also  in  general  uncertain.  Despite  the 
uncertainty  in  the  location  of  the  equilibria,  the  objective  in 
terms  of  control  design  centers  around  stabilization  of  some 
equilibrium  condition.  Typically,  one  expands  /(•,•)  about 
the  operating  point  of  interest,  say  xa,  and  then  applies 
linear  feedback  design  techniques  to  the  linearized  model. 
Static  state  feedback,  however,  does  not  apply  to  problems 
in  which  the  dynamics  and  the  targeted  operating  point 
are  uncertain.  Moreover,  static  state  feedback  changes  the 
operating  conditions  of  the  open-loop  system.  This  results 


in  wasted  control  effort  and  may  also  result  in  degrading 
system  performance. 

To  overcome  these  problems,  washout  filters  have  been 
used  in  many  applications  (e.g.,  [1],  [2],  [3],  [4],  [5],  [6], 
[7],  [8]).  A  washout  filter  (also  sometimes  called  a  washout 
circuit)  is  a  high  pass  filter  that  washes  out  (rejects)  steady 
state  inputs,  while  passing  transient  inputs  [1],  The  main 
benefit  of  using  washout  filters  is  that  all  the  equilibrium 
points  of  the  open-loop  system  are  preserved  (i.e.,  their 
location  isn’t  changed).  Thus,  one  can  concentrate  on  the 
design  of  controllers  emphasizing  the  increase  in  perfor¬ 
mance  achieved  for  a  particular  operating  point,  without 
the  potential  for  affecting  the  location  of  other  equilibria.  In 
addition,  washout  filters  facilitate  automatic  following  of  a 
targeted  operating  point,  which  results  in  vanishing  control 
energy  once  stabilization  is  achieved  and  steady  state  is 
reached. 

Although  washout  filters  have  been  successfully  used  in 
many  control  applications,  there  is  no  systematic  way  for 
choosing  the  constants  of  the  washout  filters  and  the  control 
parameters.  Recently,  Bazanella,  Kokotovic  and  Silva  [9] 
proposed  a  technique  to  control  continuous-time  systems 
with  unknown  operating  point.  The  operating  point  (or 
equilibrium  point)  was  treated  as  an  uncertain  parameter  and 
a  certainty  equivalence  adaptive  controller  was  proposed.  In 
this  work,  we  discuss  benefits  and  limitations  of  washout 
filter-aided  feedback  for  both  continuous-time  and  discrete¬ 
time  systems  .  We  also  discuss  extensions  of  washout  filter- 
aided  feedback  to  overcome  the  limitations  of  washout 
filters  and  at  the  same  time  maintain  their  benefits.  Our 
extensions  are  similar  to  that  of  [9],  although  we  do  not 
invoke  a  singular  perturbation  framework. 

The  paper  proceeds  as  follows.  In  Sec.  II,  we  discuss 
washout  filters  for  both  continuous-time  and  discrete-time 
systems.  In  Sec.  Ill,  we  discuss  linear  washout  filter- 
aided  feedback  control  and  present  limitations  of  feedback 
through  stable  washout  filters.  In  Sec.  IV,  we  discuss 
delayed  feedback  control  for  discrete-time  systems  and  its 
relation  to  washout  filter-aided  feedback.  In  Secs.  V  and  VI, 
generalizations  of  washout  filters  are  presented. 


II.  Washout  Filters 


A  washout  filter  is  a  high  pass  filter  that  washes  out  (re¬ 
jects)  steady  state  inputs,  while  passing  transient  inputs  [1], 
In  continuous-time  setting,  the  transfer  function  of  a  typical 
washout  filter  is 


G(s) 


y(s )  =  5 

x(s)  s  +  d 


1  - 


d 

s  +  d 


(3) 


Here,  d  is  the  reciprocal  of  the  filter  time  constant  which  is 
positive  for  a  stable  filter  and  negative  for  an  unstable  filter. 
With  the  notation 

2(s)  :=  (4) 

s  +  d 


the  dynamics  of  the  filter  can  be  written  as 


z  =  x  —  dz ,  (5) 

along  with  the  output  equation 

■y  =  x  —  dz.  (6) 


In  discrete-time,  the  dynamics  of  a  washout  filter  can  be 
written  as 


z(k  +  1)  =  x(k)  +  (1  -  d)z(k),  (7) 

along  with  the  output  equation 

y{k)  =  x(k)  —  dz{k).  (8) 

For  a  stable  washout  filter,  the  filter  constant  satisfies  0  < 

d<  2. 

Note  that  the  output  of  the  washout  filter  (for  both 
continuous-time  and  discrete-time  cases)  vanishes  in  steady 
state.  Therefore,  using  washout  filters  in  feedback  control 
does  not  move  the  equilibrium  points  of  the  open-loop 
system.  As  will  be  discussed  below,  there  are  limitations  in 
using  stable  washout  filters  in  feedback  control,  and  some  of 
these  limitations  can  be  overcome  using  unstable  washout 
filters. 

III.  Linear  Feedback  Through  Washout  Filters 
Below,  we  consider  linear  feedback  through  washout 
filters  for  both  continuous-time  and  discrete-time  systems, 
and  we  mention  limitations  of  using  stable  washout  filters. 
Some  of  these  limitations,  such  as  Lemma  3,  are  being 
reported  in  the  current  literature,  although  the  results  date 
back  to  the  thesis  of  H.-C.  Lee  [4],  The  results  for  the 
discrete-time  case  are  new. 

A.  Continuous-time  case 

Suppose  x0  is  an  unstable  operating  condition  for  sys¬ 
tem  (1).  In  a  small  neighborhood  of  xa,  system  (1)  can  be 
rewritten  as 


x  =  Ax  +  bu+  h(x,  u)  (9) 


evaluated  at  x0 ,  b  is  the  derivative  of  /  with  respect  to  u 
evaluated  at  xQ,  and  //.(•,  •)  represents  higher  order  terms, 
i.e.,  ft.(0,0)  =  0  and  =  q. 

Next,  washout  filters  are  used  in  the  feedback  loop.  The 
dynamic  equations  of  the  washout  filters  can  be  written  as 

n 

Zi  dj Zj:  +  ^  '  CijXj  (10) 

3=1 

where  z,  is  the  state  of  the  zth  washout  filter,  i  =  1,  •  •  • ,  to, 
and  to  <  n  is  a  positive  integer.  Note  that  (10),  where 
more  than  one  state  is  used  as  an  input  to  the  washout 
filter,  is  more  general  than  (5).  The  relationship  between 
the  operating  point  of  interest  of  the  open-loop  system  and 
the  operating  point  of  the  washout  filters  is  as  follows: 


1 

=  ~j~  y^'CjjXpj 

1  3  =  1 


(ID 


In  vector  form,  the  closed-loop  system  can  therefore  be 
written  as 


A  0  x 
C  D  )\  z 


b 

0 


h{x ,  u) 
0 


(12) 


where  C  =  [ dj ]  is  an  to  x  n  matrix,  which  consists  of 
nonzero  row  vectors,  D  =  diag(di),  i  =  1,  •  •  • ,  to. 

The  control  input  u  is  taken  as  a  linear  function  of  the 
washout  filter’  outputs  obtained  from  the  right  side  of  (10) 


IJi  d.7  Zj  ^  ^  Cjj  Xj .  (13) 

3=1 


The  following  two  lemmas  give  general  guidelines  for 
choosing  the  matrices  C  and  D  based  on  controllability 
considerations. 


Lemma  1:  ([4])  If  any  two  diagonal  entries  of  the  matrix 
D  are  the  same,  the  linearization  of  the  closed-loop  sys¬ 
tem  (12)  is  not  controllable  regardless  of  the  controllability 
of  the  pair  (A,  b). 

Proof:  From  the  Popov-Belevitch-Hautus  (PBH)  rank 
test  [10],  the  linearization  of  system  (12)  is  controllable 
if  and  only  if 


/  XI  -A  0  b 
V  C  XI  -  D  0 


71  +  TO 


for  each  complex  number  A.  Here,  p  denotes  the  rank  of  a 
matrix.  Letting  Ai  be  an  eigenvalue  of  D  with  multiplicity 
greater  than  one,  we  have 


0 

Ai  I-D 


<  TO  —  1. 


where  x  now  denotes  x  —  x0  (is  the  state  vector  referred 
to  x0),  u  is  a  scalar  input,  A  is  the  Jacobian  matrix  of  / 


Since 


P 


Ai  I- A  b 
C  0 


<  n  + 1, 


we  have 


P 


A J-A  0  b 
C  \J-D  0 


+P 


0 

Ai7  -  D 


<  P 


Ai I- A  b  \ 
C  0  ) 


<  n  +  m. 


Thus,  the  linearization  of  the  closed-loop  system  is  not 
controllable.  ■ 


with  di  >0,  i  =  !,•••,  to.  Next, 


det(^c)  =  det  ^  o  )  (  /  0  ) 

=  det  (  —bKD  A  +  bKC  ) 

=  det(-D)  det  (A  +  bKC  -  bKDD~lC) 
=  det(— D)  det(A) 


Note  that  controllability  of  the  closed-loop  system  (12) 
does  not  imply  that  the  eigenvalues  of  system  (9)  can  be 
arbitrarily  assigned  by  feedback  through  washout  filters. 


where  the  next  to  last  equality  follows  by  the  Schur  comple¬ 
ment.  Suppose  that  A  has  an  odd  number,  say  q,  of  unstable 
eigenvalues.  Then 


Lemma  2:  ([4])  Suppose  that  Ai  is  an  eigenvalue  of  both 
A  and  /),  and  that 


/  Ai  I- A  b 
{  C  0 


<  n. 


(14) 


Then,  the  linearization  of  the  closed-loop  system  (12)  is  not 
controllable. 

Proof:  Using  the  PBH  test. 


f  XJ-A  0  b 
V  C  \J-D  0 


<  n  +  to  —  1. 


(15) 


Thus,  the  PBH  fails  and  we  conclude  that  the  closed-loop 
system  is  uncontrollable.  ■ 


Since  washout  filter-aided  feedback  can  be  viewed  as  a 
form  of  output  feedback  (see  Appendix  B),  where  the  out¬ 
puts  of  the  washout  filters  instead  of  the  open-loop  system 
states  are  used  in  the  feedback,  some  of  the  capabilities  of 
direct  state  feedback  are  lost.  This  is  due  to  the  restriction 
that  di  ^  0.  The  following  lemma  summarizes  some  of  the 
capability  limitations  of  feedback  through  stable  washout 
filters. 


Lemma  3:  ([4])  If  A  has  an  odd  number  of  eigenvalues 
with  positive  real  part,  then  (9)  cannot  be  stabilized  using 
stable  washout  filters.  This  holds  even  if  the  eigenvalues  of 
A  with  positive  real  part  are  linearly  controllable. 

Proof:  Only  linear  feedback  control  is  considered  since 
nonlinear  terms  in  the  feedback  control  would  not  change 
the  linearization  of  the  system.  Using  linear  feedback  u  = 
Ky,  where  K  is  a  1  x  to  vector  and  y  is  the  vector  of 
washout  filter  outputs,  results  in  a  closed-loop  system  with 
linearization 


where 


A  +  bKC 
C 


=:  Ar 


D  = 


(  d\  0 

0  C?2 

Vo  0 


-bKD 

-D 


S' 


dm  / 


(16) 


(17) 


s?'gn(det(Ac))  =  sign(det(—D))sign(det(A)) 

=  (-l)m(-l)n-« 

=  (-l)n+m(-l)« 

By  way  of  contradiction,  suppose  that  Ac  has  no  unstable 
eigenvalues.  Then,  sign(det(Ac))  =  (— l)ra+m,  which 
contradicts  (18).  Thus,  the  closed-loop  system  possesses  at 
least  one  unstable  eigenvalue  and  cannot  be  stabilized  using 
stable  washout  filters.  ■ 

Lemma  3  implies  that  if  the  linearization  of  the  open-loop 
system  possesses  an  odd  number  of  unstable  eigenvalues, 
then  in  order  to  stabilize  the  system,  it  is  necessary  to  use 
an  odd  number  of  unstable  washout  filters  in  the  feedback 
loop. 

Corollary  1:  If  the  open-loop  system  possesses  a  zero 
eigenvalue,  it  cannot  be  moved  using  washout  filter-aided 
feedback. 

Proof:  Follows  from  the  proof  of  Lemma  3.  ■ 

B.  Discrete-time  case 

Suppose  x0  is  an  unstable  operating  condition  for  sys¬ 
tem  (2).  In  a  small  neighborhood  of  xa ,  system  (2)  can  be 
rewritten  as 


x(k  +  1)  =  Ax(k)  +  bu(k)  +  h(x(k),u(k))  (18) 


where  x  now  denotes  x  —  x0  (is  the  state  vector  referred 
to  x0),  u  is  a  scalar  input,  A  is  the  Jacobian  matrix  of  / 
evaluated  at  xa,  b  is  the  derivative  of  f  with  respect  to  u 
evaluated  at  x0,  and  represents  higher  order  terms, 

i.e.,  ft.(0,0)  =  0  and  =  0. 

Next,  washout  filters  are  used  in  the  feedback  loop.  The 
dynamic  equations  of  the  washout  filters  can  be  written  as 

n 

Zi(k  +  1)  =  (1  -  dj)zj(k)  +  yy  CjjXj(k)  (19) 

i= i 

where  Zi  is  the  state  of  the  zth  washout  filter,  i  =  1,  •  •  • ,  to, 
and  m  <  n  is  a  positive  integer.  The  relationship  between 
the  operating  point  of  the  open-loop  system  and  the  oper¬ 
ating  point  of  the  washout  filters  is  as  follows: 


- 


j  cijX0j 

J=1 


(20) 


In  vector  form,  the  closed-loop  system  can  therefore  be 
written  as 


(  x(k  + 1)  \  _  /  A  0  \  f  x(k) 

\  z{k  +  1)  )~  VC  I-D)\  z(k ) 
+  ^  h(x(k),u(k))  ^ 


where  C  =  [c-ij]  is  an  to  x  n  matrix,  which  consists  of 
nonzero  row  vectors,  D  =  diag(di),  i  =  1,  •  •  • ,  to. 

The  control  input  u  is  taken  as  a  linear  function  of  the 
washout  filter’  outputs  obtained  from  the  right  side  of  (19) 


Vi(k)  =  —djZj(k)  +  y^Cfj-a ;j(k).  (22) 

3  = 1 

The  following  two  lemmas  give  general  guidelines  for 
choosing  the  matrices  C  and  D  based  on  controllability 
considerations.  The  results  are  analogous  to  the  continuous¬ 
time  results  presented  in  the  previous  section. 


Lemma  4:  If  any  two  diagonal  entries  of  the  matrix  D  are 
the  same,  the  linearization  of  the  closed-loop  system  (21)  is 
not  controllable  regardless  of  the  controllability  of  the  pair 

{AM- 

Proof:  From  the  PBH  rank  test,  the  linearization  of  sys¬ 
tem  (21)  is  controllable  if  and  only  if 

/  XI  -  A  0  b  \ 

p{  C  XI-(I-D)  0  )=n  +  m 


for  each  complex  number  A.  Here,  p  denotes  the  rank  of 
a  matrix.  Let  di  be  an  eigenvalue  of  D  with  multiplicity 
greater  than  one.  Letting  Ai  =  1  —  d\,  we  have 


Since 


XJ-il-D) 


Ai  I- A  b 
C  0 


<  m  —  1. 


<n  +  1, 


we  have 

/  Ai  I- A  0  b\  (  Ail-  A  b\ 

V  C  XJ  -  {I  -  D)  0  )  ~  P\  C  0  ) 


P 


M  X\I  —  (I  —  D)  <n  +  m- 


Thus,  the  linearization  of  the  closed-loop  system  is  not 
controllable.  ■ 


Note  that  controllability  of  the  closed-loop  system  (21) 
does  not  imply  that  the  eigenvalues  of  system  (18)  can  be 
arbitrarily  assigned  by  feedback  through  washout  filters. 


Lemma  5:  Suppose  that  Ai  is  an  eigenvalue  of  both  A 
and  I  —  D,  and  that 


/  AiJ-A  b 
p\  C  0 


<  n. 


(23) 


Then,  the  linearization  of  the  closed-loop  system  (21)  is  not 
controllable. 


Proof:  Using  the  PBH  test. 


/  Ai  I- A  0  b 

V  C  Ai/  —  (I  —  D)  0 


<  n  +  m  —  1.  (24) 


Thus,  the  PBH  fails  and  we  conclude  that  the  closed-loop 
system  is  uncontrollable.  ■ 


Since  washout  filter-aided  feedback  can  be  viewed  as 
a  form  of  output  feedback  (see  Appendix  B),  some  of 
the  capabilities  of  direct  state  feedback  are  lost.  This  is 
due  to  the  restriction  that  di  /  0.  The  following  lemma 
summarizes  some  of  the  capability  limitations  of  feedback 
through  stable  washout  filters. 

Lemma  6:  If  A  possesses  an  odd  number  of  real  eigen¬ 
values  (counting  multiplicities)  in  (l,oo)  (i.e.,  if  det(J  — 
A)  <  0)  then  it  cannot  be  stabilized  using  stable  washout 
filters. 

Proof:  Only  linear  feedback  control  is  considered  since 
nonlinear  terms  in  the  feedback  control  would  not  change 
the  linearization  of  the  system.  Using  linear  feedback 
u(k)  =  gy(k),  where  g  is  a  1  x  m  vector  and  y(k)  is 
the  vector  of  washout  filter  outputs,  results  in  a  closed-loop 
system  with  linearization 


(  x{k+l)  \ 

V  z{k+l)  ) 


where 


D 


A  +  bgC  —bgD  \ 
C  I-D  ) 


(  di  0  •••  0  \ 

I  0  d2  ■■■  0 


\  0  0  •  •  •  dm  ) 


with  di  £  (0,2),  «  =  !,•••,  to.  Next, 


x(k)  \ 
z(k)  ) 


det  (I  —  Ac)  =  det 
=  det 


0  I  \  ,T  .  J  0  I 
I  0  )  A°}\  I  0 

D  -C 

bgD  I  —  A  —  bgC 

=  det(Z?)det(/  —  A  —  bgC  +  bgDD~1C) 
=  det (D)  det (/  —  A) 


where  the  next  to  last  equality  follows  by  the  Schur  com¬ 
plement.  If  all  the  washout  filters  are  stable,  i.e.,  di  £  (0, 2), 
then  det(D)  =  d±d2  . . .  dn  >0.  Thus,  det (/  —  Ac)  <  0  if 
det(/  —  A)  <  0.  Hence,  the  closed-loop  system  possesses 
an  odd  number  of  real  unstable  eigenvalues  in  (l,oo).  ■ 

Lemma  6  implies  that  if  the  open-loop  state  dynamics 
matrix  A  possesses  an  odd  number  of  real  eigenvalues  in 
(1,  oo),  then  in  order  to  stabilize  the  system,  an  odd  number 
of  unstable  washout  filters  (with  di  <  0)  must  be  used  in 
the  feedback  loop. 

Corollary  2:  If  the  linearization  of  the  open-loop  state 
dynamics  matrix  A  possesses  an  eigenvalue  of  1  (i.e.,  I  —  A 


is  singular),  then  this  eigenvalue  cannot  be  moved  using 
washout  filter-aided  feedback. 

Proof:  Follows  from  the  proof  of  Lemma  6.  ■ 

IV.  Delayed  Feedback  Control  as  a  Special 
Case  of  Washout  Filters 

Delayed  feedback  control  (DFC)  was  introduced  by  Pyra- 
gas  [11]  as  a  technique  for  control  of  chaos.  Since  its 
introduction,  DFC  has  been  used  in  many  applications.  It 
has  been  shown  in  [12]  that  DFC  for  discrete-time  systems 
has  two  limitations.  The  first  limitation  is  known  as  the  odd 
number  of  real  eigenvalues  greater  than  1  limitation.  That  is, 
DFC  cannot  be  used  to  stabilize  systems  whose  linearization 
possesses  an  odd  number  of  real  eigenvalues  in  (l,oo). 
The  second  limitation  is  that  DFC  can  be  used  to  stabilize 
only  a  class  of  unstable  systems;  it  cannot  stabilize  highly 
unstable  systems  [13].  For  example,  for  the  one  dimensional 
map  (25)  with  system  dynamics  coefficient  a,  DFC  can 
be  used  to  stabilize  the  map  if  and  only  if  —3  <  a  < 
— 1  (see  below).  A  discussion  of  two-dimensional  discrete 
time  systems  that  can  be  stabilized  using  DFC  is  given 
in  [12].  Many  researchers  proposed  different  techniques  to 
overcome  the  odd  number  limitation  (e.g.,[14],  [15],  [16], 
[17]).  Extended  delayed  feedback  control  (EDFC),  where 
many  previous  states  of  the  system  are  used  in  the  feedback, 
was  proposed  to  extend  the  range  of  systems  that  can  be 
stabilized  (e.g.,  [13]).  However,  the  analysis  of  the  EDFC 
method  is  cumbersome  and  it  also  suffers  from  the  odd 
number  limitation  as  DFC  [18]. 

In  this  work,  we  note  that  delayed  feedback  control 
for  discrete-time  systems  is  a  special  case  of  washout 
filter-aided  feedback  (it  corresponds  to  washout  filter-aided 
feedback  with  all  washout  filters’  constants  equal  1).  We 
then  proceed  to  show  that  the  limitations  of  DFC  mentioned 
above  can  be  overcome  by  the  use  of  washout  filter-aided 
feedback,  giving  rise  to  the  possibility  of  stabilizing  a  much 
larger  class  of  systems  than  is  possible  with  DFC. 

To  illustrate  DFC  and  its  limitations,  consider  the  simple 
one  dimensional  discrete-time  system 

x(k  +  1)  =  ax(k)  +  u(k)  (25) 


J  is  Schur  stable  if  and  only  if 

-1  <  7  <  1 
0+7  <  1+7 
0  +  7  >  —1  —  7 

The  second  inequality  implies  that  a  <  1  and  the  first  and 
third  inequalities  imply  that  a  >  —3.  Of  course,  the  system 
is  already  stable  if  —  1  <  a  <  1.  Thus,  DFC  can  stabilize 
an  unstable  plant  (25)  iff  —3  <  a  <  —1. 

Next,  we  show  that  DFC  for  the  one  dimensional  sys¬ 
tem  (25)  is  a  special  case  of  washout  filter-aided  feedback. 
We  then  show  that  using  washout  filter-aided  feedback,  any 
one-dimensional  system  can  be  stabilized. 

Consider  the  same  system  as  before  but  with  washout 
filter-aided  feedback: 


x(k  + 1) 

=  ax{k )  +  u(k) 

(28) 

z(k  + 1) 

=  x(k)  +  (1  —  d)z{k) 

(29) 

u(k ) 

=  7  (x(k)  —  dz(k)) 

(30) 

Here  d  is  the  washout  filter  constant  and  7  is  the  control 
gain.  Note  that  the  washout  filter-aided  feedback  reduces  to 
DFC  by  setting  d  =  1. 

The  fixed  point  of  the  closed-loop  system  is  asymptoti¬ 
cally  stable  if  the  eigenvalues  of  the  closed-loop  system  are 
within  the  unit  circle.  Note  that  if  \a\  <  1,  the  uncontrolled 
system  (i.e.,  (28)  with  u(k )  =  0)  is  asymptotically  stable. 

Proposition  1:  Denote  the  system  dynamics  coefficient 
by  a.  Suppose  that  a  <  —1  or  a  >  1,  i.e.,  the  open-loop 
system  is  unstable. 

Case  1:  If  a  <  —1,  then  a  stabilizing  washout  filter- 
aided  feedback  exists  (using  a  stable  washout  filter).  In¬ 
deed,  washout  filter-aided  linear  feedbacks  with  gain  7  and 
washout  filter  constant  d  satisfying 

—  1  —  a  +  ^(1  +  a)  <  7  <  1  —  a(l  —  d),  0  <  d  <  —  (31) 

are  stabilizing. 

Case  2:  If  a  >  1,  then  a  stabilizing  washout  filter- 
aided  feedback  exists  (using  an  unstable  washout  filter). 
Indeed,  washout  filter-aided  linear  feedbacks  with  gain  7 
and  washout  filter  constant  d  satisfying 


where  u(k)  is  the  control  input.  In  DFC,  the  control  is  taken 
to  be  u(k)  =  7 (x(k)  —  x(k  —  1)),  where  7  is  the  control 
gain. 

The  closed-loop  system  can  be  written  as 


—  1  —  a  +  ^(1  +  a)  <  7  <  1  —  a(  1  —  d),  ^  <  d  <  0  (32) 

are  stabilizing. 

Proof:  See  Appendix  A. 


x(k  +  1)  =  ax(k)  +  7  (x(k)  —  z(k))  (26) 

z{k  +  1)  =  x{k)  (27) 


A  necessary  and  sufficient  condition  for  DFC  to  be  sta¬ 
bilizing  is  —3  <  a  <  1  [19],  [12],  [20].  This  can  be 

seen  from  (26)-(27)  as  follows:  The  linearization  is  J  := 

a  +  7  —7 


1 


0 


.  By  Jury’s  test  for  second  order  systems. 


For  the  n-dimensional  case,  DFC  takes  the  form 

x(k  +  1)  =  Ax{x)  +  bu{k)  (33) 

u{k)  =  g{x{k)  —  x{k  —  1))  (34) 

This  corresponds  to  washout  filter-aided  feedback  with  n 
uncoupled  washout  filters  (i.e.,  in  Eq.  (19),  =  1  if  i  = 

j  and  Cij  =  0  if  i  ^  j  )  and  all  washout  filters  having 


constants  equal  to  1  (i.e.,  in  Eq.  (19),  di  =  1,  i  =!,■■■  ,ri): 


x(k  +  1) 

=  Ax(x)  +  bu(k) 

z(k  +  1) 

=  x(k) 

y{k) 

=  x(k)  —  z(k) 

u(k ) 

=  g(x(k)  -  z{k)) 

Although  the  case  of  one-dimensional  systems  was  han¬ 
dled  in  a  systematic  way  above,  systems  of  higher  dimen¬ 
sion  do  not  lend  themselves  to  a  general  analytical  design 
procedure.  After  presenting  a  two-dimensional  example,  we 
will  proceed  in  the  next  section  to  give  a  generalization  of 
washout  filter-aided  feedback  that  does  permit  development 
of  a  systematic  design  method. 


Example  1:  Consider  the  two-dimensional  map  [21] 

(  x\(k  +  1)  \  _  /  1.9  1  \  /  xi(k)  \  _  f  x\(k)  \ 
V  x2(k+l)  )  \  0.5  0  j  \  x2{k)  /  \  o  y 

+  (  J  )  u(k)  (35) 


The  uncontrolled  system  (35)  has  three  fixed  points  x0±  = 

^ ,  xo2  =  ^  )  and  X°3  =  ~x°2’  and  indeed 

displays  chaotic  motion  (see  [21]).  The  fixed  point  x0i  is 
unstable:  the  eigenvalues  of  the  linearization  at  the  origin 
are  Ai  =  2.1343  and  A2  =  —0.2343.  Since  A2  >  1,  the  ori¬ 
gin  cannot  be  stabilized  using  DFC  nor  using  stable  washout 
filters.  We  will  show  that  the  origin  can  be  stabilized  using 
one  unstable  washout  filter 


The  authors  of  [21]  use  dynamic  delayed  feedback  con¬ 
trol  and  reduced  order  dynamic  DFC  to  stabilize  the  origin 
of  (35).  Dynamic  DFC  results  in  tripling  the  dimension  of 
the  system  and  reduced  order  dynamic  DFC  results  in  a 
closed-loop  system  dimension  more  than  twice  the  dimen¬ 
sion  of  the  open-loop  system.  On  the  other  hand,  washout 
filter-aided  feedback  results  in  increasing  the  dimension  to 
at  most  twice  the  dimension  of  the  open-loop  system.  In 
the  example  above,  only  one  washout  filter  was  needed  to 
stabilize  the  system. 


k 


Fig.  1.  Time  series  (with  initial  condition  (1.4725,0.0557))  of  (a)  x\  (b) 
X2  and  (c)  control  input  u.  The  control  is  applied  when  the  trajectory  of 
the  open-loop  system  enters  the  neighborhood  {x  =  (x\,X2)  E  • 
||(r||  <  0.1}  of  the  origin. 


z(k+  1)  =  xi(k)  +  (1  —  d)z(k)  (36) 

y(k)  =  Xi(k)  —  dz(k)  (37) 

where  z(k)  is  the  washout  filter  state,  y(k )  is  the  washout 
filter  output  and  d  is  the  washout  filter  constant  the  value 
of  which  determines  the  stability  of  the  washout  filter.  Fet 
u(k)  =  'yy(k),  where  7  is  the  control  gain  to  be  chosen. 

The  closed-loop  system  takes  the  form 

/  xx(k  +  1)  \  /  1-9  +  7  1  -7 d  \  /  xi(k)  \ 

x2(k+l)  =  0.5  0  d  x2(k) 

\  z(k  +  1)/  \  1  0  1  —  d  J  \  z(k)  J 

x\{k)  \ 

0  (38) 

0  J 

The  origin  is  locally  asymptotically  stable  if  the  lineariza¬ 
tion  of  (38)  is  asymptotically  stable.  Since  the  uncontrolled 
system  has  one  eigenvalue  greater  than  1,  we  need  to  choose 
d  such  that  the  washout  filter  is  unstable  but  the  closed- 
loop  system  is  stable.  It  is  an  easy  calculation  to  show 
that  a  controller  with  d  =  —0.05  and  7  =  —1.8  stabilizes 
the  origin  locally.  With  such  parameters,  the  eigenvalues 
of  the  Jacobian  matrix  of  the  closed-loop  system  are 
{—0.6345, 0.8923  ±  jO. 1767}.  Figure  1  demonstrates  the 
effectiveness  of  the  controller. 


V.  Generalization  of  Continuous-Time  Washout 
Filter-Aided  Feedback 

Next,  we  consider  a  generalization  of  washout  filters  in 
which  the  individual  washout  filters  are  coupled  through 
a  constant  coupling  matrix.  Consider  system  (1)  with  x0 
as  the  operating  condition.  System  (1)  can  be  rewritten  as 
follows  in  a  small  neighborhood  of  xa: 

x  =  Ax  +  Bu  +  h(x,  it)  (39) 

We  are  interested  in  designing  a  control  law  that  stabilizes 
this  system  while  maintaining  all  its  equilibrium  points. 

The  generalized  washout  filter-aided  feedback  proposed 
here  results  in  the  closed-loop  system 


Ax  +  Bu  +  h(x,  u) 

(40) 

P{x  -  z) 

(41) 

K  (x  —  z) 

(42) 

Here  P  is  a  nonsingular  matrix  and  K  is  a  feedback  gain 
matrix.  Since  in  steady  state,  the  control  input  vanishes  (i.e., 
u  =  0),  the  equilibrium  points  of  the  open-loop  system  are 
not  shifted  by  this  type  of  feedback  control.  Suppose  that 
the  pair  (A,  B)  is  stabilizable.  Are  there  matrices  K  and  P 
such  that  the  closed-loop  system  is  stable? 


To  answer  this  question,  we  consider  the  effect  of  matri¬ 
ces  K  and  P  on  the  linearization  of  the  closed-loop  system 


x 

z 


(  A  +  BK 

-BI<  \ 

(  x  \ 

{  P 

P  ) 

z 

is  Hurwitz.  Recall  that  eigenvalues  are  preserved  under 
similarity  transformations. 


Let  /  — 


P- 


.  Then  we  have 


Aci  TiAcTl  1 


A  +  BK  - BKP  \ 
I  -P  ) 


(47) 


Proposition  2:  ([9])  The  determinant  of  the  closed-loop 
state  dynamics  matrix  Ac  satisfies 

det(Ac)  =  det(A)  det(— P)  (44) 

Proof:  Using  the  fact  that  similarity  transformations  do  not 
change  the  eigenvalues  and  the  Schur  complement  of  a 
matrix,  we  have  that 

det(Ac)  =  det  ^  o  )  (  /  0  ) 

=  det  (  —BK  A  +  BK  ) 

=  det(-P)  det  (A  +  BK  -  BKPP"1) 

=  det(— P)  det(A) 


Corollary  3:  If  the  matrix  A  has  a  zero  eigenvalue,  then 
the  closed-loop  system  state  dynamics  matrix  Ac  will  also 
have  a  zero  eigenvalue. 

Proof:  Follows  from  Proposition  2.  ■ 

Since  this  type  of  feedback  doesn’t  shift  equilibria,  it 
shouldn’t  be  surprising  that  it  can’t  modify  a  zero  eigenvalue 
(this  would  also  modify  any  stationary  bifurcation  in  the 
system). 

The  following  result  gives  some  conditions  on  the  con¬ 
troller  matrix  P  for  the  controller  to  be  stabilizing.  This 
result  is  akin  to  Lemma  3  pertaining  to  washout  filter-aided 
feedback.  In  words,  the  result  means  that  if  the  open-loop 
system  possesses  an  odd  number  of  unstable  eigenvalues, 
then  a  necessary  condition  for  the  closed-loop  system  to  be 
stable  is  that  the  controller  must  also  have  an  odd  number 
of  unstable  eigenvalues. 

Lemma  7:  [9]  Let  the  number  of  unstable  eigenvalues 
of  A  be  odd.  Then,  for  the  closed-loop  matrix  Ac  to  be 
Hurwitz,  —  P  must  also  have  an  odd  number  of  unstable 
eigenvalues. 

Proof:  Recall  that  det(Ac)  =  det  (A)  det(— P).  Thus, 

sign(det(Ac))  =  sign(det(A))sign(det(—P))  (45) 

Let  q  be  the  number  of  unstable  eigenvalue  of  A,  and  r  be 
the  number  of  unstable  eigenvalues  of  —P.  Then  using  the 
fact  that  Ac  is  Hurwitz,  Eq.  (45)  becomes 


Next,  let  T2  = 


I  M 
0  I 


.  It  is  easy  to  see  that  T2  1  = 


I  -M 
0  I 


.  Applying  the  transformation  T2  to  Ac\  gives 


-4C2  :=  T2AciT2 

(a  +  BK  +  M  -AM- BKM-M2-BKP-MP 
I  -M  -  P 


Consider  the  (1,2)  block  term  of  Ac2.  Suppose  P  =  eP\ 
and  M  =  Mo  +  eMi  +  0(e2)  with  e  >  0  and  sufficiently 
small.  It  is  straightforward  to  show,  after  setting  the  (1,2) 
block  of  Ac 2  to  zero,  i.e.. 


AM  +  BI\M  +  M2  +  BKP  +  MP  =  0  (48) 


and  collecting  terms  with  same  power  in  e,  that  0(1)  terms: 

{A  +  BK  +  M0)  M0  =  0.  (49) 

This  holds  if  M0  =  0  or  M0  =  -A-  BK.  Taking  M0  = 
—A  —  BK  and  finding  the  e1  terms  gives 

Mi  {A  +  BK)  +  AP1  =  0  (50) 

Since  A  +  BK  can  be  guaranteed  invertible  (by  restricting 
K  so  that  0  ^  a(A  +  BK)),  we  find  that  Mi  =  —AP\(A  + 
BK)"1.  Since  Mi  can  be  determined  uniquely  through 
matrix  inversion,  it  is  clear  that  the  Implicit  Function 
Theorem  implies  that  (48)  has  a  locally  unique  solution 
M(e)  =  Mo+eMi+0(e2)  near  Mq.  Therefore,  M  =  Mo+ 
eMi  +  0(e2)  =  -A  -  BK  -  eAPi (A  +  BK)"1  +  0(e2). 
Substituting  M  and  P  in  Ac2  yields 

_  f  -eAP^A  +  BK)"1 +  0(e2)  0  \ 

c2~V  1  Ac2(2,2)  ) 

where  Ac2(2,2)  =  A  +  BK  +  e(AP1{A  +  BK)"1  -  Px)  + 
0(e2). 

Assume  that  A  has  no  zero  eigenvalues.  To  make  Ac2 
Hurwitz,  we  need  to  choose  Pi  such  that  —APi(A+BK)~1 
is  Hurwitz.  Clearly  such  a  Pi  exists  (e.g..  Pi  =  ,4  ~ 1  ( .4  + 
BK)).  Also  we  need  to  choose  K  such  that  A  +  BK  is 
Hurwitz  with  eigenvalues  away  from  zero  such  that  the 
perturbation  e(APi(A+  BK)"1  —  Pi)  does  not  cause  the 
eigenvalues  of  A  +  BI\  to  become  unstable.  Such  a  K  is 
guaranteed  to  exist  since  the  pair  (A,B)  is  assumed  to  be 
stabilizable. 


(-1)2™  =  (-l)"-‘?(-l)n-’~  =  (— l)2n(— 1)9+,’(46) 

Thus,  q  +  r  must  be  even.  ■ 

We  will  show  that  if  A  is  nonsingular  and  ( A ,  B) 
is  stabilizable,  then  there  is  a  pair  P,  K  such  that  Ac 


Proposition  3:  Consider  the  closed-loop  system  (43). 
Suppose  that  the  matrix  A  has  no  eigenvalues  at  0.  Suppose 
also  that  the  pair  (A,  B)  is  stabilizable.  Then  there  exists 
a  P  €  Rnxn  and  K  g  Rmxn  sucj1  ^xJ,x^)T  is 
asymptotically  stable  equilibrium  point  of  (43). 


VI.  Generalization  of  Discrete-Time  Washout 
Filter-Aided  Feedback 

The  results  of  this  section  are  counterparts  of  the 
continuous-time  results  of  the  previous  section  for  the 
discrete-time  case.  Consider  system  (2)  with  x0  as  the 
operating  condition.  System  (2)  can  be  rewritten  as  follows 
in  a  small  neighborhood  of  x0: 

x(k  +  1)  =  Ax(k)  +  Bu(k)  +  h(x(k),  u(k))  (51) 

We  are  interested  in  designing  a  control  law  that  stabilizes 
this  system  while  maintaining  all  its  equilibrium  points.  The 
generalized  washout  filter-aided  feedback  proposed  here 
results  in  the  closed-loop  system 

x(k  +1)  =  Ax{k)  +  Bu{k )  +  h(x(k),u(k))  (52) 

z(k  +  1)  =  Px(k)  +  (/  -  P)z{k)  (53) 

u(k)  =  G(x(k)  -  z(k ))  (54) 

Here  P  is  a  nonsingular  matrix  and  G  is  a  feedback  gain 
matrix.  Since  in  steady  state,  the  control  input  vanishes  (i.e., 
u  =  0),  the  equilibrium  points  of  the  open-loop  system  are 
not  shifted  by  this  type  of  feedback  control.  Suppose  that 
the  pair  ( A ,  B)  is  stabilizable.  Are  there  matrices  G  and  P 
such  that  the  closed-loop  system  is  stable? 

To  answer  this  question,  we  consider  the  effect  of  matri¬ 
ces  P  and  G  on  the  linearization  of  the  closed-loop  system 

(  x{k  +  l)\  _  f  A  +  BG  -BG  \  (  x(k)  \ 

\z(k+l)J  ^  P  I-P  )\  z{k)  ) 

-=*($>) 

Proposition  4:  The  determinant  of  I  —  Ac  satisfies 

det  (7  —  Ac)  =  det(7  —  A)  det(P)  (56) 
Proof:  Using  the  fact  that  similarity  transformations  do  not 
change  the  eigenvalues  and  the  Schur  complement  of  a 
matrix,  we  have  that 

det  (7  —  Ac )  = 


Corollary  4:  If  the  open-loop  matrix  A  has  an  eigenvalue 
of  1  (i.e.,  if  (I  —  A)  is  singular),  this  eigenvalue  cannot  be 
shifted  using  this  type  of  dynamic  feedback. 

Proof:  Follows  from  Proposition  4. 

Lemma  8:  Let  the  number  of  unstable  eigenvalues  of  A 
that  are  real  and  greater  than  1  be  odd  (i.e.,  det(7— A)  <  0). 
Then,  for  the  closed-loop  state  dynamics  matrix  Ac  to  be 
Schur  stable,  I  —  P  must  also  have  an  odd  number  of  real 
eigenvalues  greater  than  1  in  value. 


det  (  /  o  )  det('J  ~  Ac ^  det  (  /  o  ) 

det  (  BG  I  —  A  —  BG  ) 
det(P)  det  (7  -A-BG+  BGP^P) 
det(P)  det(7  —  A) 


Proof:  We  show  that  for  Ac  to  be  Schur  stable,  it  is  neces¬ 
sary  that  the  parity  of  the  number  of  real  eigenvalues  of  P 
that  are  negative  (equivalently,  number  of  real  eigenvalues 
of  (7  —  P)  that  are  greater  than  1)  be  equal  to  that  of  A. 
Recall  that  det(7  —  Ac )  =  det(7  —  A)  det(P).  Thus, 

sign(det(I  —  Ac ))  =  sign(det(I  —  A))sign(det(P))  (57) 

Let  q  be  the  number  of  real  eigenvalue  of  A  with  value 
greater  than  1,  and  r  be  the  number  of  negative  real 
eigenvalues  of  P.  Then  using  the  fact  that  Ac  is  Schur 
stable,  Eq.  (57)  becomes 

1  =  (— 1)9(— l)r  =  (-l)9+r  (58) 

Thus,  q  +  r  must  be  even.  ■ 

We  will  show  that  there  exist  a  P,  G  such  that  Ac  is 
Schur  stable.  Recall  that  eigenvalues  are  preserved  under 
similarity  transformations. 

Let  Ti  =  ^  d  .  Then  we  have 

Ac  i  '■=  TiAcTi  1 

A  +  BG  - BGP  \ 

I  I -P  )  ' 

Next,  let  T2  =  o  ^  = 

^ .  Applying  the  transformation  P2  to  Ac\  gives 

Aq2  '■=  T2AciT2  1 

_  (  A  +  BG  +  M  Ac2(1,  2)  \ 

^  7  -A7  +  7-P  ) 

where 

Ac2(1,  2)  =— AM — BGM — M2  —  BGP + M — M P.  (59) 

Consider  the  block  term  ,4,:2  (1.2).  Suppose  P  =  ePi  and 
M  =  +  eMi  +  0(e2)  with  e  >  0  and  sufficiently  small. 

It  is  straightforward  to  show,  after  setting  the  Ac2(l,2)  to 
zero  and  collecting  terms  with  same  power  in  e  that  the 
0(1)  terms  yield 

(A  +  BG  -  7  +  M0)M0  =  0  (60) 

This  holds  if  Mq  =  0  or  Mq  =  —A  —  BG  +  I.  Taking 
Mo  =  —A  —  BG  +  I  and  finding  the  e1  terms  gives 

M1(A  + BG  -  I)  +  (A- I)P1  =0  (61) 

Since  A  +  BG  —  I  can  be  guaranteed  nonsingular  (by 
restricting  I\  so  that  1  ^  cr(A  +  BG)),  we  find  that 
Mi  =  —{A  —  I)Pi{A  +  BG  —  7)-1.  Since  Mi  can  be 
determined  uniquely  through  matrix  inversion,  it  is  clear 
that  the  Implicit  Function  Theorem  implies  that  (59)  has  a 
locally  unique  solution  M(e)  =  Mq  +  eMi  +  0(e2)  near 
M0.  Therefore,  M  =  M0  +  eMi  +  0(e2)  =  —A  —  BG  + 
I  —  e(A  —  I)Pi{A  +  BG  —  7)-1  +  0(e2).  Substituting  M 
and  P  in  Ac2  yields 

I  —  e(A  —  I)Pi(A  +  BG-  I)-1  +  0{e2) 

I 


f  I  -M 
l  0  7 


0 

j4c(2,  2) 


where 


2 


Ac{  2,2)  =  A  +  BG  +  eHA-^P^A  +  BG-I)-1  -Pi) 

+0(e2). 

Assume  that  I  —  A  is  nonsingular  (i.e.,  1  ^  cr(A)).  To 
make  AC2  Schur  stable,  we  need  to  choose  Pi  such  that 
I  e(A  —  I)Pi(A+BG— I)”1  is  Schur  stable.  Clearly  such 
a  P±  exists.  Also  we  need  to  choose  G  such  that  A  +  BG 
is  Schur  stable  with  eigenvalues  away  from  1  such  that  the 
perturbation  e((A  —  I) Pi  (A  +  BG  —  J)-1  —  Pi)  does  not 
cause  the  eigenvalues  of  A  +  BG  to  become  unstable.  Such 
a  G  is  guaranteed  to  exist  since  the  pair  (A,  B)  is  assumed 
to  be  stabilizable. 

Proposition  5:  Consider  the  closed-loop  system  (55). 
Suppose  that  the  matrix  I  —  A  is  nonsingular.  Suppose 
also  that  the  pair  ( A ,  B)  is  stabilizable.  Then  there  exists  a 
P  G  Rnxn,  a  G  G  Rmxn  and  an  e  >  0  such  that  Ve  G  (0,  e], 
(xJ,x^)T  is  an  asymptotically  stable  fixed  point  of  (55). 

Example  2:  (Design  of  a  stabilizing  controller  using  the 
generalized  washout  filter  calculations) 

We  revisit  Example  1  using  the  generalized  washout  filter 
design  calculations  above.  We  choose  the  gain  vector  G  so 
that  A+bG  is  Schur  stable.  A  stabilizing  control  gain  vector 
is  G  =  [-1.6343  -  0.7657].  Choosing 

Pi  =  (A  —  I)~1(A  +  BG  —  I) 
f  -0.1674  -0.5469  \ 

V  -°-5837  0.7265  ) 

and  e  =  0.1  yields  P  =  ePi,  and 


( 

A  +  BG 

—BG  \ 

\ 

P 

d-P) 

/ 

0.2657 

0.2343 

1.6343 

0.7657  \ 

0.5000 

0 

0 

0 

-0.0167 

-0.0547 

1.0167 

0.0547 

V 

-0.0584 

0.0727 

0.0584 

0.9273  / 

The  eigenvalues  of  Ac  are 

{-0.2343,  0.7277,  0.8164,  0.9000}.  Thus,  the  closed- 
loop  system  is  asymptotically  stable.  Figure  2  demonstrates 
the  effectiveness  of  the  controller.  Note  that  the  control 
input  vanishes  after  stabilization  of  the  origin  is  achieved. 

Appendix  A  Proof  of  Proposition  1 
The  Jacobian  of  the  closed-loop  system  (28)-(30)  is 
t (  a  +  7  -7 d  \ 

V  1  1  -d) 

Let  r  :=trace( J)  =  a  +  1  —  d  +  +,  <5  :=det(J)  =  a(l  — 
d)  +  7.  The  fixed  point  of  the  closed-loop  system  (28)-(30) 
is  asymptotically  stable  if  both  eigenvalues  of  J  are  within 
the  unit  circle. 

The  characteristic  equation  of  J  is  given  by 
p( A)  :=  A2  -  (a  +  7  +  1  —  d) A  +  a(l  —  d)  +  7 


Fig.  2.  Time  series  (with  initial  condition  (0.3, -0.6))  of  (a)  x\  (b)  X2  and 
(c)  control  input  u.  The  control  is  applied  when  the  trajectory  of  the  open- 
loop  system  enters  the  neighborhood  { x  =  (xi,X2)  E  :  \\x\\  <  0.15} 
of  the  origin. 


By  the  Jury’s  test  for  second  order  systems,  both  eigenval¬ 


ues  are  within  the  unit  circle  if  and  only  if 

-1  <  p(0)  <  1  (62) 

p(  1)  >  0  (63) 

p(-l)  >  0  (64) 

Conditions  (62)-(64)  are  equivalent  to 

—  1  <  a(l  —  d)  +  7  <  1  (65) 

d{l-a)  >  0  (66) 

2  +  2a  +  27  —  d  —  ad  >  0  (67) 


respectively. 

Case  1:  a  <  —  1 

Let  d  G  (0,2).  This  corresponds  to  a  stable  washout  filter. 
Then,  inequality  (66)  is  trivially  satisfied  since  d  is  positive 
and  a  <  —1  by  hypothesis.  Inequalities  (65)  and  (67) 
translate  to  an  explicit  condition  on  7  and  d  as  follows: 

max  |— 1  —  a(l  —  d),  —  1  —  a+^(l  +  a)|<7< 

1  -  a(l  -  d)  (68) 

Now,  max  {— 1  —  a(l  —  d),  —  1  —  a  +  |(1  +  a)}  = 

— 1  —  a  +  |(1  +  a),  which  is  seen  as  follows: 

—1  —  a(l  —  d)  <  —  1  —  a  +  —  (1  +  a) 

,  d.,  , 

ad  <  -(1  +  a) 

a  <  1  (tme  by  hypothesis;  a  <-l) 

Hence,  inequality  (68)  reduces  to: 

—  1  —  a  +  ^(1  +  a)  <  7  <  1  —  a(l  —  d)  (69) 

Finally,  for  a  7  satisfying  (69)  to  exist,  the  upper  limit 
in  (69)  must  be  greater  than  the  lower  limit.  This  is  shown 


to  be  true  as  follows: 
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1  —  a(l  —  d)  >  —  1  —  a  +  -(1  +  a) 

d 

2  -f-  ad  >  —  (1  ^0 

^  ad  a 

2  +  t>  7, 

d  < 


1  —  a 


Case  2:  a  >  1 

Similar  to  the  previous  case.  From  inequality  (66),  d(  1  — 
a)  >  0  if  and  only  if  d  <  0.  Inequalities  (65)  and  (67) 
translate  to  an  explicit  condition  on  7  and  d  as  follows: 


max  —  1  —  a(l  —  d),  —  1  — a+-(l  +  a)j<7< 

1  -  a(l  -  d)  (70) 


which  reduces  to 

—  1  —  a  +  ^(1  +  a)  <  7  <  1  —  a(l  —  d)  (71) 

For  a  7  to  exist,  the  upper  limit  in  (71)  must  be  greater  than 
the  lower  limit.  This  implies  that  d>  This  completes 
the  proof.  ■ 


Appendix  B  Washout  Filter-Aided  Feedback  as 
an  Output  Feedback 

Below,  we  show  that  washout  filter-aided  feedback  can 
be  viewed  as  an  output  feedback.  Consider 

x(k  +  1)  =  Ax{k)  +  bu(k)  (72) 

with  washout  filter-aided  feedback 

z(k  +  1)  =  Cx{k)  +  (I  —  D)z(k)  (73) 

y(k)  =  Cx{k)  —  Dz{k)  (74) 

Here,  C  =  [c^]  is  an  m  x  n  matrix  and 
D  =diag{di,  c?2,  •  •  • ,  dm},  m  <  n.  Let  the  control 
input  be 

u(k)  =  gy(k)  =  g{Cx{k)  -  Dz(k))  (75) 

where  the  row  vector  g  is  the  gain  vector. 

The  closed-loop  system  (72)-(75)  can  be  written  as 

-  (s:i!) 

=  (c  i-d){:u )+(!)■« 

=:  Ax(k)  +  bu(k ) 

»(*>  =  <c  -D>  (*<*)) 

=  :  Cx{k) 
u(k)  =  gy(k) 
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